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Abstract 

The weh-known supersymmetric constructions such as Witten's super- 
symmetric quantum mechanics, Spiridonov-Rubakov parasupersymmetric 
quantum mechanics, and higher-derivative SUSY of Andrianov et al. are 
extended to the nonstationary Schrodinger equation. All these construc- 
tions are based on the time-dependent Darboux transformation. The su- 
peralgebra over the conventional Lie algebra is constructed. Examples of 
time-dependent exactly solvable potentials are given. 

1. The supersymmetric quantum mechanics is in a permanent intensive devel- 
opment since the Witten papers Q. One can cite the iV-extended supersymmetric 
quantum mechanics [|[, parasupersymmetric quantum mechanics [^, and higher 
order derivative supersymmetric quantum mechanics The field of supersym- 
metric quantum mechanics is recently reviewed in We want to point out that 
all above mentioned constructions are valid for the time independent Hamilto- 
nians and if one restrict oneself by the stationary solutions of the Schrodinger 
equation. Hence, these constructions can be referred to the stationary super- 
symmetric quantum mechanics and the nonstationary one needs to be developed. 
We hope that this report gives a stimulus for the further developments in this 
domain. 

2. The nonstationary supersymmetric quantum mechanics is based on the 
nonstationary Darboux transformation just in the same way as the stationary 
one ^ is based on the conventional Darboux transformation |^]. 

Let us consider two time-dependent Schrodinger equations 

{idt-Ho)ip{x,t) = o, dt = d/dt, Ho = -d^,-Vo{x,t), dl = dA, (i) 

{idt- Hi)ip{x,t) Hi = -dl-Vi{x,t), x e R, t G (2) 

Here — Vo(x,i) is a potential energy and R = [a,b] is the interval for variable x 
which can be both finite and infinite. If the Schrodinger operators for Eqs. (1) 
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and (2) are related by the intertwining relation 



L{idt - Ho) = {idt - Hi)L, (3) 

where L is a linear operator, called transformation operator, the functions ip and 
ip are related as follows: f = Lip if Lip ^ 0. 

If (idt—Ho) and {idt— Hi) are self-adjoint (in the sense of some scalar product) 
the equation (3) implies 

L+{idt-Hi) = {idt-Ho)L+, (4) 

where the superscript plus sign (+) is used to denote the operator adjoint to L, 
and Eqs. (1) and (2) become "peer". It follows from Eqs. (3) and (4) that 
So = L~^L commutes with (idt — Hq) and si = LL~^ commutes with (idt — Hi) 
and consequently sq is a symmetry operator for the initial equation (1) and si is 
a symmetry operator for the final one (2). 

The constructions such as in Eq. (3) are well-known in mathematics and are 
intensively investigated since the Delsart's paper |^. The most significant results 
obtained with the help of the transformation operators concern the backscattering 



problem in quantum mechanics |1C] and its application for the solving of the 
nonlinear equations [pi]]. 

3. We now assume operator L to be a differential of the first degree in dx 
with smooth coefficients depending on both variables x and t. We should not 
include in L the derivative dt since it, being found from equation (1), transforms 
L into the second-order differential operator. In this case the operator L and the 
real potential difference A{x,t) = Vi{x,t) — Vo{x,t) are completely defined by a 
function u{x,t) called transformation function 

L = Li{-Ux/u + dx), (5) 
Li = Li{t) =exp{2jdtlin(logu)xx), (6) 
A = {\og\u\'^)xx- (7) 

To obtain a real potential difference we should impose the reality condition ^] 
on the function u 

(lognK),,, = 0, (8) 

where the asterisk implies the complex conjugation. 

In the majority of cases of physical interest we can introduce the Hilbert 
space structure Lq{R) in the space of solutions of the equation (1) with the scalar 
product appropriately defined. Symmetry operator sq = L^L, being symmetric, 
may be extended up to self-adjoint in the appropriate Hilbert space and can have 
either discrete spectrum or continuous one. Since Lu = (see Eq. (5)), the 
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function u is an eigen function of this operator corresponding to zero eigen value. 
Hence, u is one of the eigen functions of operator L+L = h — a. In general, h is 
a self-adjoint integral of motion for the initial Schrodinger equation (1) which in 
particular case (if Vq does not depend on t) may be equal to the Hamiltonian Hq. 

With the help of the transformation operator L, just in the same way that 
in the conventional supersymmetric quantum mechanics [||, 0], we can construct 
the supercharge operators 



Q 




L'ip{x, t) 

ger equations (1) and (2) may now be rewritten in supersymmetric form 



acting on two-component wavefunctions 

Ly now be rewritten 

{ildt-n)-^{x,t) = 0, 



(9) 



Two Schrodin- 



(10) 



where / is 2 x 2 identity matrix and 7i 



is a super hamiltonian. 



' Ho 

Since so = L^L and si = are symmetry operators for equations (1) and (2) 

respectively, the superoperator S = ^ LL+ ^ symmetry operator 

for the equation (10). The operators Q, Q^, and S form a well-known superalge- 
bra [H, 0. There is a single difference, namely, instead of the Hamiltonians we see 
other integrals of motion of the equation (1) in its construction. When h is the 
initial Hamiltonian above constructions coincide with known ones. This is the 
reason to call the transformation (5), (6), time- dependent Darboux transformation 

i- 

4. With the help of the other eigen functions of operator h we can perform 
the chain of Darboux transformations and construct the parasuperalgebra in full 
analogy with papers If in this chain we eliminate all intermediate operators 
h and express the final operator L in terms of particular solutions Ui of the ini- 
tial equation (1), we obtain higher-derivative nonstationary quantum mechanics 
similar to the stationary one [S]. In this case 



L = L^^^ = LN{t)W-\ui, . . . . 



UN 



Ulx U2x 



1 

dx 



u 



(N) „(N) 



Ix 



2x 



d: 



N 



(11) 



where W stands for the conventional symbol of the Wronskian of the functions 
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Ul, . . . ,un. For the real function Lfyf(t) we have 

Liv(t)=exp{2 j dtlm[logW{ui,...,UN)]xx}- (12) 
The reahty condition (8) takes now the form 

For the potential difference we obtain 

AN = {log\W{uu...,UN)\'),,. (14) 

We can recognise in formula (11) the generalisation of the known Crum-Krein 
formula |12, ^] to the nonstationary case. Note that the condition (13) is more 



feeble than the reality condition (8) imposed on every function Ui. Thus, we can 
construct the higher-derivative supersymmetry with the self-adjoint final Hamil- 
tonian even if the intermediate Hamiltonians are not self-adjoint (so-called ir- 
reducible case described for stationary equation in Ref. 4). The basic relation 
for the time-dependent polynomial supersymmetric quantum mechanics is the 
following factorisation properties 

N N 

L+L = l[iho-a), LL+ = l[ihi-Ci) (15) 

1=1 1=1 

first obtained for the stationary case in Ref. 4. The Cj in Eqs. (15) are eigen 
values corresponding to the eigen functions Ui of the operator ho. 

5. To obtain a regular potential difference by the Darboux transformation 
(5)-(7) the transformation function u should be nodeless. In the space Lq{R) 
a single nodeless eigen function of the operator h exists (if it has a discrete 
spectrum). This function is the ground state function. Beyond the space Lq{R) 
there are many nodeless eigen functions of the operator h suitable for use as 
the transformation functions. They should have eigen values a < Eq {eq being 
the lowest eigenvalue of h). In this case the discrete spectra of the symmetry 
operators h = L^L + a and h = LL~^ + a differ by a single level and we have 
unbroken supersymmetry Q. Every bounded state of the superoperator S, except 
for its ground state, is double degenerate. 

We will now describe unexpected peculiarities in the breakdown of the super- 
symmetry in the higher-derivative supersymmetric quantum mechanics. These 
peculiarities (as far as we know) are not discussed in the available literature. 

The Darboux transformation (5,6) being performed with the discrete spec- 
trum function Unix, t) of the operator h [it has n zeros in the interval (a, b)] gives 
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the potential difference with n poles. Solutions obtained with the help of the 
transformation operator (5) are not square integrable functions in [a,b]. Never- 
theless, the second transformation with the transformation function n„_|_i(x,t), 
having n + 1 zeros in (a, b) removes all singularities and the transformation op- 
erator of the second degree L^^^ = L2L1, where L12 are the first order Darboux 
transformation operators, is well defined. This fact reflects the known property of 
the Wronskian constructed from the functions li^. belonging to the space L'^{R) 



[13|: the Wronskian W{uk^ , ■ ■ ■ , ^fciv) conserves its sign if for all k = 0, 1, 2, . . . 
the inequality {k — ki){k — /C2) ■ ■ ■ {k — k^) > holds. In particular, the functions 
Mfc. may be two by two juxtaposed functions. The discrete eigenvalues a^. of the 
operator h are absent in the spectrum of transformed operator H. This signifies 
that the ground state level of the superoperator S is two-fold degenerate and the 
excited states constructed with the help of the functions Uk^ are nondegenerate. 
Furthermore, these states are annihilated by the operators Q and whereas the 
ground state is annihilated only by the one of these operators. It should be noted 
that this property remains valid for the stationary case, i.e., in the conventional 
supersymmetric quantum mechanics. 

6. Differential symmetry operators for the stationary Schrodinger equation 
are Hamiltonian and its polynomial functions. Symmetry algebra of the nonsta- 
tionary Schrodinger equation is more rich than for the stationary case. We can 
use this algebra to construct a superalgebra. For this purpose we should define 
an operator inverse to L. 

The equation (5) implies Lu = 0. Choose the transformation function u such 
that the absolute value of u~^{x, t) is square-integrable in the interval R and the 
condition (8) is satisfied. Then for every ip £ Lq(R) we have if = Lip £ Ll(R), 
but the set L^i^R) = {ip : ip = Lip, tp £ Lg(i?)} does not span the whole space 
L\{R). The function (y9o(x,t) = [Li{t)u* {x,t)\~^ € L\{R) can not be obtain 
by the action of the operator L on any ip G Lq{R). If we designate by L\q{R) the 
linear hull of the function lpq then L\{R) = LIq{R) © Lf^(i?). 

Choose as the transformation function the function v = lpq and define the 
following integral operator acting from Lf^(i?) to Lq{R) 

Mip{x,t) = [Li{t)v*{x,t)]-^ r v*{y,t)ip{y,t)dy. (16) 



The straightforward calculation persuades that LMip = ip for all ip G L\^{R) and 
the condition v G L\q{R) implies MLip = ip for all ip G Lq{R) operator M, hence, 
IS inverse to L : M = L-^ and we have one-to-one correspondence between the 
spaces Lq{R) and Lii{R). 

If in the space Lq{R) the symmetry operators Qi forming n-dimensional Lie 
algebra G with the structural constants flj : [gi,gj] = fljQi are defined and this 
space is invariant under the action of these operators then in the space L\y{R) 
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we can define the operators gi = LgiM and this space will be invariant under the 
action of all gi. Furthermore these operators form a basis for the n-dimensional 
Lie algebra G with the same structural constants flj. Let Tq be the space of 
two-component wave functions with the basis ^+{x,t) = ijj{x,t)e+ and 

^'_(x,i) = Lip{x,t)e-, V £ L'q{R), and e+ = (J), e_ = (5*). In the space Tq we 
can define operators 

which form a Lie algebra isomorphic to G. The operators Gi are symmetry 
operators for the supersymmetric equation (10). Besides the operators Gi in 
the space Tq the following operators can be defined: Pq = La~ , Qi = giL^^a^ 

where ~ 0^''^^~^0 0^' '^^^^^ operators are evidently nilpotent: 

Pq = 0, = 0, and {Qi, Qj} = QiQj + QjQi = 0. Furthermore, we can find by 
the direct calculations that {Po,Qi} = Gi and the generalised Jacobi identities 
are fulfilled. The operators Gi,Qi,Po, i = 1,2, ...,n, hence, form a basis for 
2n + 1-dimensional Lie superalgebra sG. 

Note that since h = L^L + a £ G we have for the operator S introduced in 
sec. 3, S G sG and Q, £ sG. 

7. Examples. Consider first the simplest case of a free particle: vo{x,t) = 0. 
Choose the following solutions of the initial Schrodinger equation (1) [p!4| : 

7px{x,t) = {l + t')-^/^exp[ix't/{4: + At') + iXarctant]Qxiz), (18) 

Z = x(l +t2)-V2^ 

where Qxiz) is the parabolic cylinder functions satisfying the equation Q'x{z) — 
{z'^ /A+X)Qx{z) = with A being an arbitrary parameter (a separation constant). 
For A = n + 1/2, n = 0,1,2,..., the functions Qx{z) are expressed via the 
Hermite polynomials Qn+i/2{^) = exp(z^/4)i/„(i2;/\/2). The reality condition 
(8) is satisfied for all real A. Functions (18) being nodeless for A = n + 1/2 and 
for even n are suitable for use as transformation functions. Formula (7) gives the 
new potential 

)(x,t) = -(1 + + 4ki2k - 1)^?^ - 8k'(^^^Y), 

^ q2k{z) ^ q2k{z) ' ' 

where q^^z) = {-ifHek{iz), Hek{z) = 2~^l'Hk{z/ ^). 

The same functions for odd n are nodeless in the semiaxis (0, cxd) and with 
their help we obtain the following time-dependent exactly solvable potential 

(2fc+l) _ ,^ , .2^-l(. , .r/r,, , iN92fc-l(^ , -1^21^ <l2k{z) ^^2^ 



-{l + t')-Ul + ^k{2k + lf-^l^^-2{2k + lf( 

^ q2k+i{z) ^q2k+i{z) 



:Vl + t' e (0,oo). 
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The functions (18) for A = —n — 1/2 form a discrete basis set in Lq{M.^). The 
double Daxboux transformation with juxtaposed functions tx„ = V'-n-i/2 and 
Un+i gives a regular potential of the form 

4«'^^+i)(.x,t) = -2{l + t^)mz)/Mz) - {J'^{z)IUz)f - 1], 

Jk{z) = T{k + 1) E r-i(s + l)Hel{z) = kJk-i{z) + Heliz), 

s=0 

Jo{z) = 1, Ji{z) = z^ + l, J2{z) = + 3, . . . . 

These are the potentials which correspond to the supersymmetric model with 
two-fold degenerate cigcnvahie of the supcropcrator S except for the eigenvahics 
corresponding to the functions u„ and Un+i- When n > the ground state of 5* 
is two-fold degenerate and nondegenerate eigenvalues are situated in the middle 
of the spectrum of S 

It is not difficult to establish that the Wronskian constructed from two func- 
tions (18) with Ai = m + 1/2 and A2 = / + 1/2 for m = 0, 2,4, . . . and / = 
m + l,m + 3, . . . is nodeless and therefore these functions are suitable for double 
Darboux transformation. This gives the following exactly solvable potential 

vt''\x,t) = -2{l + er\l + d^logfU^)/dz^), 
fml{z) = qmiz)qi+i{z) - qi{z)qra+i{z). 

We will cite one example for harmonic oscillator potential as well: Hq = 
-dl + iJ'x^, Hoi^n = (2ra + l)*n, i^n = Hn{^/iox) exp ( - iuj{2n + l)t - 

n = 0,1,2, If we choose the following nonstationary solution of the initial 

Schrodinger equation (1) as the transformation function: 

u{x, t) = sin~-'-/^(2a;t) cos h{\x/ sin 2ut) x 
exp[i(a;a;^ - \^/uj) cot(2a;t)/2] ^ Lo(K^), A G M^ 

we obtain the nonstationary anharmonic potential of the form: vi{x, t) = w^x^ — 
2A2 sin- V2 (2a;i )sech2 {\x / sin 2ujt)). 

First author was supported by the Russian Foundation for Fundamental Re- 
search. 

References 

[1] E. Witten, Nucl. Phys., B 188 (1981) 513; ibid., B 202 (1982) 253. 

[2] A.I. Pashnev, Sov. J. Theor. Math. Phys., 60 (1986) 311; V.P. Berezovoy 
and A.I. Pashnev, Sov. J. Theor. Math. Phys., 78 (1989) 289; V.P. Berezovoy 
and A.I. Pashnev, Z. Phys. 51 (1991) 525. 



7 



[3] V.A. Rubakov and V.P. Spiridonov, Mod. Phys. Lett. A3 (1988) 1337; A.A. 
Andrianov and M.V. loffe, Phys. Lett. B 255 (1991) 543; J. Beckers, N. 
Debergh, and A.G. Nikitin, Mod. Phys. Lett. 7 (1992) 1609. 

[4] A.A. Andrianov, M.V. lofFe, and V.P. Spiridonov, Phys. Lett. A 174 (1993) 
273; A.A. Andrianov, M.V. loffe, and D.N. Nishnianidze, Polynomial SUSY 
in Quantum Mechanics and Second Derivative Darhoux Transformation, 
preprint SPbU-IP-94-05 (1994); A.A. Andrianov, F. Canata, J.-P. Dedonder, 
and M.V. loffe. Second Order Derivative Supersymmetry and Scattering 
Problem, preprint SPbU-IP-94-03 (1994). 

[5] F. Cooper, A. Khare, and U. Sukhtame, Phys. Rep. 251 (1995) 267. 

[6] V.G. Bagrov, B.F. Samsonov, and L.A. Shekojan, Izv. Vyssh. Uchebn. 
Zaved. Fizika No 7 (1995) 59; V.G. Bagrov and B.F. Samsonov, Phys. Lett. 
A (in press). 

[7] A.A. Andrianov, N.V. Borisov, M.V. loffe, and I.M. Fides, Sov. J. Thcor. 
Math. Phys. 61 (1984) 17; C.V. Sukumar, J. Phys. A 18 (1985) 2917; ibid., 
2937; V.G. Bagrov and B.F. Samsonov, Sov. J. Theor. Math. Phys. 104 
(1995) 356. 

[8] G. Darboux, Compt. Rend. Acad. Sci. 94 (1882) 1456; G. Darboux, Legons 
sur la theorie generale des surfaces et les application geometrique du calcul 
infinitesimal. Deuxieme partie, Paris, Gauthier-Villars et fils (1889). 

[9] J. Delsart, J. Math. Pures et Appl. 17 (1938) 213; J. Delsart and J.L. Lions, 
Comment. Math. Helv. 32 (1957) 113. 

[10] Z.S. Agranovich and V.A. Marchenko, Backscattering Problem, Kharkov 
(1969); B.M. Levitan, Inverse Sturm- Liouville Problems, Nauka, Moscow 
(1984); L.D. Faddeev, Usp. Mat. Nauk 105 (1959) 57. 

[11] V.B. Matveev and M.A. Salle, Darboux transformations and solitons, Berlin, 
Springer (1991); F. Calogero and A. Degasperis, Spectral Transform and 
Solitons, Amsterdam-New York-Oxford; R.K. Dodd, J.C. Eilbeck, J.D. Gib- 
bon, and H.C. Morris, Solitons and Nonlinear Wave Equations, Academic 
Press; V.E. Zakharov, Backscattering Method, Solitons. Ed. by R. Bulla and 
P.M. Kodri, Moscow (1983) 270. 

[12] M. Crum, Quart. J. Math. 6 (1955) 263. 

[13] M.G. Krein, Dokl. Akad. Nauk SSSR 113 (1957) 970. 

[14] W. Miller Jr., Symmetry and Separation of Variables. Massachusetts (1977). 



8 



